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; ABSTRACT 

, We explore empirical relations between three different populations of high-redshift 

^ ' galaxies and their hosting dark halos employing the halo model approach. Specifically 

I we consider LBGs (Lyman-break galaxies) at z ~ 4 and at z ~ 5, and LAEs (Lyman- 

■ Alpha emitters) at z ~ 4.86, all from the Subaru Deep Field survey extending over an 

' area of about 600 arcmin^. We adopt a halo occupation function (HOF) prescription to 

parameterize the properties of their hosting halos and the efficiency of halo-dependent 
star formation. We find that the two LBG samples are well described by the halo 
model with an appropriate HOF. Comparing the model predictions with the observed 
number densities and the angular correlation functions for those galaxies, we obtain 
• constraints on properties of their hosting halos. A typical mass of hosting halos for 

Qh' LBGs is 5 X l0^^h~^MQ and the expected number of LBGs per halo is ~ 0.5, therefore 

Q . there is an approximate one-to-one correspondence between halos and LBGs. We also 

^ • find a sign of the minimum mass of LBG hosting halos decreasing with time, although 

[ its statistical significance is not strong. We discuss implications of these findings on the 

i star formation history of LBGs. On the other hand, for LAEs, our simple HOF pre- 

^ ' scription fails to reproduce simultaneously the observed angular correlation function 

. 1-H , and the number density. In particular, a very high amplitude of the correlation func- 

' tion on scales larger than 120 arcsec cannot be easily reconciled by the HOF model; 

[ a set of parameters which account for this high correlation amplitude on large scales 

predict either excessive clustering on small scales or a much smaller number density 
than observed. While this difficulty might imply either that the distribution of LAEs 
within hosting halos differs from that of dark matter, or that the strong large-scale 
correlation is due to the existence of an unusual, large overdense region, and so the 
survey region is not a representative of the z ~ 5 universe, the definite answer should 
wait for a much wider survey of LAEs at high redshifts. 
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galaxies: formation — dark matter 



o 
cn 
o 



1 INTRODUCTION 

Multi-band color selection techniques (Steidel et al. 1996; 
1998; Madau et al. 1996; Cowie & Hu 1998; Ouchi et al. 2001, 
2003a) have significantly increased high-redshift galaxy cat- 
alogs both in quality and in size. Since those high-2 galaxies 
are naturally expected to be progenitors of the present-day 
galaxies, their statistical analysis is of fundamental impor- 



tance in understanding the formation and evolution history 
of galaxies. Actually, recent large high-z galaxy catalogs al- 
low one to estimate their luminosity functions and spatial 
correlation functions at difi'erent z with a reasonable accu- 
racy (e.g., Steidel et al. 1999; Adelberger et al. 1998; Gi- 
avalisco & Dickinson 2001; Ouchi et al. 2001, 2003a; Porciani 
& Giavalisco 2002). 

In the standard scenario of structure formation, it is 
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thought that dark matter halos are first formed via the grav- 
itational amplification of initial small density fluctuations. 
Subsequently baryonic gas trapped in a gravitational po- 
tential of the dark matter halo becomes sufficiently dense 
to cool and to form stars, and such initially small systems 
experience repeated mergers to form larger galaxies. The 
formation process of halos is determined by gravity only, 
thus it is well understood from TV-body simulations and also 
from simple but relevant analytical approximations such as 
the Press-Schechter model (Press & Schechter, 1974) and 
its extensions (Bond et al. 1991; Bower 1991; Lacy & Cole 
1993; Mo & White 1996; Sheth, Mo & Tormen 2001; Sheth 
& Tormen 2002). The formation of galaxies, on the other 
hand, involves many complicated processes including hy- 
drodynamics, radiative processes, star formation, and su- 
pernova feedback, and thus it is hard to solve in a reliable 
manner even using state-of-the-art numerical simulations. 
Therefore, a simplified or empirical model that describes 
essential physical processes is still a valuable tool to un- 
derstand basic elements of the formation and evolution of 
galaxies. 

This is why we attempt in what follows to apply an em- 
pirical parameterized model that relates the galaxy number 
distribution to mass of the hosting dark matter halo. To be 
more specific, we explore statistical relations between two 
populations of high-z galaxies and their hosting dark mat- 
ter halos: (i) Lyman break galaxies (hereafter LBGs) which 
are isolated in a color-color diagram due to their UV contin- 
uum depression (Steidel et al. 1996; 1998; Madau et al. 1996; 
Adelberger et al. 1998; Ouchi et al. 2001) and (ii) Lyman 
Q emitters (LAEs) which are identified due to their strong 
Lyman ct emission from narrow band imaging (Cowie & Hu 
1998; Hu, Cowie & McMahon 1998). We consider three cat- 
alogs generated from the Subaru Deep Field survey data 
(Ouchi et al. 2001, 2003a; Shimasaku et al. 2003); LBGs at 
z ~ 4, LBGs at z ~ 5 and LAEs at z ~ 4.68. 

The major purpose of our current analysis is twofold; 
the first is to clarify the difference of hosting halos for 
LBGs and LAEs. It is well known that these two popula- 
tions exhibit different statistical properties, including the 
fact that LAEs are in general fainter and smaller, and are 
more strongly clustered than LBGs (Ouchi et al. 2003a). 
These differences should retain information of their forma- 
tion processes as well as environmental effects. The second 
is to examine the difference between LBGs located at differ- 
ent redshifts (z ~ 4 and ~ 5). Combined with the previous 
analysis of LBGs at z ~ 3 by Moustakas & Somerville (2001) 
and Bullock et al. (2002), our results would provide better 
understanding of the evolution of LBGs. 

For those purposes, we employ the halo approach that 
attempts to model the spatial distribution of galaxies in a 
parameterized fashion. The key quantity in this approach is 
the halo occupation function (HOF) that describes statisti- 
cal relations between galaxies and their hosting halos. We 
adopt a simple form for HOF motivated by the semi-analytic 
galaxy formation models (Benson et al. 2000; Kauffmann et 
al. 1999) as well as by hydrodynamic simulations (Yoshikawa 
et al. 2001; White, Hernquist & Springel 2002). Combined 
with models of the halo mass function and spatial cluster- 
ing of halos, for which very accurate fitting functions are 
obtained from A^-body simulations (Jing 1998; Sheth & Tor- 
men 1999; Hamana et al. 2001; Jenkins et al. 2001), the halo 



approach predicts the spatial clustering of galaxies as well 
as their number density (Mo & Fukugita 1996; Mo, Mao & 
White 1999). Comparing these predictions with the observed 
values, we obtain constraints on the relation between galax- 
ies and their hosting halo mass for different populations of 
galaxies. The latter methodology was first attempted by Jing 
& Suto (1998) for LBGs at z ~ 3 using their halo catalogs 
from A'^-body simulations. 

The outline of this paper is as follows. Section 2 de- 
scribes models and basic equations. Section 3 summarizes 
observational data that are used to put constraints on the 
model parameters. In section 4, our results of LBGs at z ~ 4 
and ~ 5 are presented and are compared with previous re- 
sults of z ~ 3 LBGs (Moustakas & Somerville 2001; Bul- 
lock et al. 2002). We show results of LAEz5, and discuss 
their implications in section 5. Finally, section 6 is devoted 
to summary and discussion. In Appendix, we illustrate how 
the HOF parameters depend on the shape of the two-point 
correlation function. 

Throughout this paper, we adopt a flat ACDM (Cold 
Dark Matter) cosmology with the matter density Q,^ — 0.3, 
the cosmological constant Qa ~ 0.7, the Hubble constant 
Ho = 100/ikm/s/Mpc with h — 0.7, and the normalization 
of the matter power spectrum ag = 0.9. We adopt the fit- 
ting function of the CDM power spectrum of Bardeen et 
al. (1986). 



2 HALO APPROACH FOR GALAXY 
CLUSTERING 

The basic idea behind the halo model that we adopt be- 
low has a long history (Neyman & Scott 1952; Limber 1953; 
Peebles 1974, 1980; McClelland & Silk 1977; and other re- 
cent papers referred to in this section). The model predic- 
tions have been significantly improved with the recent ac- 
curate models for the mass function, the biasing and the 
density profile of dark matter halos (Seljak 2000; Peacock 
& Smith 2000; Ma & Fry 2000). This approach has been 
applied to various problems in cosmological nonlinear clus- 
tering, galaxy clustering and weak lensing correlation (e.g., 
Sheth & Jain 1997; Jing, Mo & Borner 1998; Komatsu & Ki- 
tayama 1999; Cooray, Hu & Miralda-Escude 2000; Cooray 
& Sheth 2002; Scoccimarro et al. 2001; Shu, Mao & Mo 
2001; Hamana, Yoshida & Suto 2002; Berlind & Weinberg 
2002; Bullock, Wechsler & Somerville 2002; Moustakas & 
Somerville 2002; Takada & Jain 2003; Takada & Hamana 
2003). In this section, we summarize several expressions 
which are most relevant to the current analysis. In particu- 
lar, we mainly follow the modeling of Berlind & Weinberg 
(2002), Bullock et al. (2002) and Moustakas & Somerville 
(2002) in which readers may find further details. 

We adopt a simple parametric form for the average 
number of a given galaxy population as a function of the 
hosting halo mass: 



TVJM) 



f (M/Mi) 
1° 



(M > Mn 

{ M < M„ 



(1) 



The above statistical and empirical relation is the essential 
ingredient in the current modeling characterized by the min- 
imum mass of halos which host the population of galaxies 
(Mmin), a normalization parameter which can be interpreted 
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as the critical mass above which halos typically host more 
than one galaxy (Mi; note that Mi may exceed Mmin since 
the above relation represents the statistically expected value 
of the number of galaxies), and the power-law index of the 
mass dependence of the efficiency of galaxy formation (a). 
We will put constraints on the three parameters from the 
observed number density and clustering amplitude for each 
galaxy population. In short, the number density of galaxies 
is most sensitive to Mi which changes the average number of 
galaxies per halo. The clustering amplitude on large angu- 
lar scales (> 1') is determined by the hosting halos and thus 
very sensitive to the mass of those halos, Mmin. The cluster- 
ing on smaller scales, on the other hand, depends on those 
three parameters in a fairly complicated fashion; roughly 
speaking, Mmin changes the amplitude, a, and to a lesser 
extent Mi as well, changes the slope (see Appendix). Fur- 
ther detailed discussion may be found in Berlind & Weinberg 
(2002), Bullock et al. (2002) and Moustakas & Somerville 
(2002). 

With the above relation, the number density of the cor- 
responding galaxy population at redshift z is given by 



.(^) 



dM nhaio(M,2) Ng{M), 



(2) 



where nhaio(M) denotes the halo mass function for which 
we adopt the fitting function of Sheth & Tormen (1999). 

The galaxy two-point correlation function on small 
scales is dominated by contributions of galaxy pairs lo- 
cated in the same halo. We adopt the following model 
(Bullock et al. 2002) for the mean number of galaxy pairs 
(JVg(Afg - 1))(M) within a halo of mass M: 



Nl{M) 



if iVg(M) > 1 



(iVg(7Vg-l))(M) 



Nl{M) log(4iVg(M))/log(4) 







if 1 > iVg(M) > 0.25 
otherwise. 



(3) 



The above empirical model is motivated by previous results 
from the semi-analytic galaxy formation models (Benson et 
al. 2000; Kauffmann et al. 1999) which indicate that for 
N^{M) > 1 the scatter around the mean number of galax- 
ies is Poissonian, while for N^(M) < 1 it becomes sub- 
Poissonian. 

In the framework of the halo model, the galaxy power 
spectrum consists of two contributions, one from galaxy 
pairs located in the same halo (1-halo term) and the other 
from galaxy pairs located in two different halos (2-halo 
term): 



(4) 



Assuming the linear halo bias model (Mo & White 1996), 
the 2-halo term reduces to 



Pii„(fc) 

— / dMnhaio(M)iVg(M)6(M)y(fc,M)r5) 



where Piin(fc) is the linear dark matter power spectrum, 
6(M) is the halo bias factor (we adopt the modified fit- 
ting function of Sheth & Tormen 1999), and y{k,M) is the 
Fourier transform of the halo dark matter profile normal- 
ized by its mass, y{k, M) — p{k, M)/M. See, e.g., section 3 
of Seljak (2000) for details. Here we assume that galaxies in 



halos trace the density profile of the underlying dark halos by 
Navarro, Frenk & White (1996; 1997), and adopt the mass- 
concentration parameter relation by Bullock et al. (2001) 
but with an appropriate correction (see Shimizu et al. 2003). 
Since the clustering on large scales is dominated by the 2- 
halo term, it is fairly insensitive to the assumption of galaxy 
distribution inside the hosting halo (Berlind & Weinberg 
2002). It should be noted that since j/ ~ 1 on large scales 
(e.g., scales much larger than the virial radius of halos), on 
such scales the 2-halo term can be rewritten by 

P^^ik) = bl{> M^,„)Pnn{k) (6) 

where the galaxy number weighted bias factor is defined by 

_ / dM nhaio(Af)iVg(M)b(M) 



bg{> M„ 



J dM nhaio(M)iVg(M) 



(7) 



Note that bg{> Mmin) depends on a. 
The 1-halo term is written as 



(27r)3ni. 
|y(fc,M)|*' 



dM nhaio(M) {N^iNg 



1)){M) 
(8) 

We choose p = 2 for (iVg(iVg - 1)) > 1 and p = 1 for 
(iVg(iVg - 1)) < 1 (Seljak 2000). Here we assume that in 
the limit of a small number of galaxies in one halo, one 
galaxy is located near the center of the halo. Therefore, in 
this case, the number of pairs is dominated by the central 
galaxy paired with a halo galaxy, and thus the probability 
of finding a galaxy pair is given by the single density profile 
of the galaxies within a halo. 

Once the power spectrum of the galaxy population is 
specified, one can easily compute their angular two-point 
correlation function via the Limber projection (see e.g, chap- 
ter 2 of Bartelmann & Schneider 2001): 



0(6) = Jdr q\r) j gfc Pg(fc, 



r) .h[fK{r)ek], 



(9) 



where q{r) is the normalized selection function and Jo(x) 
is the zeroth-order Bessel function of the first kind. For 
the spatially flat cosmology (fim -I- Q.a ~ 1) as we consider 
throughout the present paper, the radial function fKij) is 
equivalent to r, and r — r{z) is the radial comoving distance 
given by 

dz 



r{z) 



^0 Jo ya^TTIpTnl 



(10) 



The dependences of EOF parameters on the shape of the 
angular two-point correlation function are summarized in 
Appendix. 

For a given selection function of the observation, the 
average galaxy number density is 



dz '^^y^ q{r)ng^2{r) 
dz 



dz -^q(r) 



(11) 



where dV{r)/dz denotes the comoving volume element per 
unit solid angle: 

^^r^^z)-^- (12) 

again for the spatially flat cosmology. 
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Table 1. Summary of observed properties, here the large scale 
bias is defined by fe = ^Jw^j/wdm on scales > 100". 



Figure 1. Selection functions as a function of redshift for LBGz4 
(top), LBGz5 (middle) and LAEz5 (bottom). 
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Figure 2. Angular two-point correlation functions (upper panel) 
and the corresponding bias parameter (lower panel) for three pop- 
ulations of high-redshift galaxies, crosses, open squires, and filled 
circles represent the measured angular two-point correlation func- 
tions, •Wg{9), for LBGz4, LBGz5, and LAEz5, respectively. In the 
lower panel, points for LBGz4 and LAEz5 are slightly shifted hor- 
izontally for clarify. The measured correlation functions plotted 
have been corrected for contamination and for integration con- 
stant (see text). The model predictions for dark matter angular 
two-point correlation functions, ■w^^(8), in the ACDM cosmology 
with the same selection functions are plotted in curves; LBGz4 
(solid line), LBGz5 (dotted line) and LAEz5 (dashed line). To 
compute them, the nonlinear fitting function of the CDM power 
spectrum by Peacock & Dodds (1996) is used. The biasing pa- 
rameter is simply defined as b{9) = Wg (9)/ (9). 



Sample 


number density [h^fApc ^] 


large scale bias 


LBGz4 


(5.86 ±0.71) X 10-3 


3-4.5 


LBGz5 


(8.05 ±4.96) X 10-4 


5-7 


LAEz5 


(3.01 ± 1.94) X 10-3 


5-9 



3 DATA 

We use three different samples of galaxy populations from 
deep imaging data taken as part of the Subaru Deep Field 
(SDF) survey; LBGs at z ~ 4 (LBGz4s), LBGs at z ~ 5 
(LBGz5s) and LAEs at z ~ 4.86 (LAEz5s). Observational 
details of those samples are described in Ouchi et al. (2001, 
2003a, 2003b), and thus we briefly summarize their basic 
features relevant for our comparison. 

LBGz4s are selected from an i'-detection catalog con- 
structed from deep BRi' imaging data over a 543 arcmin^ 
area in the SDF (Ouchi et al. 2003a) . The hmiting AB mag- 
nitudes for the 3 a detection of an object in a l".8 diameter 
aperture are B = 27.8, R = 27.1 and i' — 26.9. In order to 
guarantee a reasonable level of photometric completeness, 
the j'-detection catalog is limited to i' — 26.5, correspond- 
ing to the absolute magnitude of Mnoo = —19.0-1-5 log h for 
LBGz4s. A total of 1438 LBGz4 candidates are detected in 
a range of 3.5 < z < 4.5. The selection function is shown in 
top-panel of Figure 1, which is estimated using the Monte- 
Carlo simulation on the basis of colors and redshifts of Hub- 
ble Deep Field North galaxies given in the photometric red- 
shift catalog by Furusawa et al. (2000; see Ouchi et al. 2003c, 
Ouchi 2003 for details). 

It should be noted that the selection functions plotted in 
Figure 1 are the probability distribution functions of galaxy 
redshifts in our samples. The current observational method 
does not specify the redshift of individual galaxies accurately 
except in a statistical sense. Their number density is esti- 
mated to be nLBGz4 = (5.86 ± 0.71) x IQ-'-^h^Mpc'^ . The 
angular two-point correlation function is computed by the 
procedure described in Ouchi et al. (2001; 2003d) and is 
plotted in Figure |21 

LBGz5s are selected from a ^'-detection catalog con- 
structed from deep Vi' z' imaging data over a 616 arcmin^ 
area in the SDF (Ouchi et al. 2003b). The limiting AB mag- 
nitudes are B = 27.8, V = 27.3, R = 27.1, i' = 26.9, 
z' — 26.1 for the 3 a detection in a l".8 diameter aperture. 
Again for the photometric completeness, the z'-detection 
catalog is limited to z' — 26.0, corresponding to the ab- 
solute magnitude of Afiroo = —19.7 + 51ogh, for LBGz5s. 
A total of 246 LBGz5 candidates are detected in a range 
of 4.2 < 2 < 5.2. Their number density is estimated to be 
WLBGzs = (8.05 ±4.96) X lO-^/i^Mpc-^. 

LAEz5s were first selected from the same data of 
LBGz5s but an additional observation using a narrow band- 
fiher (NB711, central wavelength of 7126 ± 4 A, FWHM 
bandwidth of 73.0 ± 0.6 A) was performed to identify LAEs 
at z ~ 4.86 (Ouchi et al. 2003a). The limiting magni- 
tude is NB711= 26.0 for the 3 a detection in a l".8 di- 
ameter aperture. The selection function of LAEz5 is ap- 
proximated by a top-hat function (bottom-panel of Fig- 
ure 1) whose shape (center and width) is determined from 
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Figure 3. Confidence contour maps derived from Ax^ for LBGz4 on tlie two-parameter plane after marginalizing over the remaining 
one parameter. Top-left panel is on Mmin-Afi, top-right on a-M\ and bottom-left on M^i^-a. A darker gray-scale indicates a lower 
Ax'^ value (thus more likely). Contour lines indicate from inner to outer Ax^ = 2.3, 6.17 and 11.8, which, if each bin of the correlation 
function is independent, correspond to 68.3, 95.4, and 99.73% confidence levels, respectively. In the present analysis, these confidence 
levels should be understood as approximate estimates. 




Figure 4. Same as Figure ISl but for LBGz5. 
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the central wavelength and FWHM of the NB711 filter 
(Ouchi et al. 2003a). A total of 87 LAEz5 candidates are 
detected over 4.83 < 2: < 4.89. Their number density is 
riLAEzs = (3.01 ± 1.94) X lO'^/i^Mpc-^ 

The angular two-point correlation function Wg{0) for 
each sample is plotted in Figure|21together with that for dark 
matter, Wdm(^), which is calculated applying the same selec- 
tion function of the observational data. The plotted correla- 
tion functions have been corrected for contamination and for 
the "integration constant" (see Groth & Peebles 1977). Pos- 
sible contaminants are mainly low-2 galaxies which match 
the selection criteria by chance, and thus we may safely as- 
sume no cross-correlation between high-z LBGs/LAEs and 
low-z contaminants and no auto-correlation in low-2 con- 
taminants. Therefore we have made contamination correc- 
tion for each sample by multiplying the observed correla- 
tion function by a factor of 1/(1 — fc)^, where the con- 
tamination rate fc is estimated to be 0.01, 0.26 and 0.4, 
for LBGz4s, LBGz5s and LAEz5s, respectively. The inte- 
gration constant is estimated to be 0.00676, 0.00637 and 
0.00675 for LBGz4, LBGz5 and LAEz5, respectively (Ouchi 
et al. 2001, 2003a, 2003b). We define the biasing parameter 
as b{6) = ^ Wg{6)/wdin{9)- For LBGs, we obtain large-scale 
(specifically, 60" < 6* < 1000") bias factors of b = 3 - 4.5 
and 5 — 7 for LBGz4 and LBGz5s, respectively, and the bias 
increases with decreasing the separation. On the other hand, 
LAEs exhibit stronger clustering on larger scales {9 > 200") 
than LBGs, while their biases on smaller scales are similar 
to those of LBGs. 



4 RESULTS ON LYMAN-BREAK GALAXIES 

4.1 Constraints on Mmin, Mi, and a 

We estimate the range of allowed values for the three pa- 
rameters, Mmin, Ml, and a by considering the following 
function constructed from the observed number density and 
the angular two-point correlation functions: 



x'(Mn,in,Mi,a) = 



[C^°'^^(gbin)-C^'"°^'"'(^bin)]' 
0-S(6'bin) 



+ 



[lognf'^ -logn^"'°'*'=']2 



(13) 



log n„ 



where (t„ and (Tiog are the statistical 1-a error in the mea- 
surements of the angular correlation function and the num- 
ber density, respectively. In the above likelihood estimator, 
we take the logarithm of the galaxy number density instead 
of the number density itself, because the predicted galaxy 
number density varies logarithmically with Mi. Note that 
although the HOF parameters can depend on time in gen- 
eral, we assume here that the three parameters are constant 
within the redshift interval of each sample for simplicity. For 
LAEz5s, this must be the case as the redshift interval is very 
small. It turns out (see t|4.3ll that the HOF parameters for 
LBGs do not change significantly over three LBG samples at 
z ~ 3, z; ~ 4 and z ~ 5. Therefore the above assumption is 
reasonable for LBGs as well. We also note that the present 
analysis does not take the cosmic variance into account al- 
though it may be important for relatively small survey vol- 
umes for those samples. Figures |3 and 2] show the map on 
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Figure 5. Comparison of the observed angular correlation func- 
tion of LBGz4s with the model prediction assuming M^i^ = 
1.6 X IO^/i-IMq, Ml = 8 X lO^^h-^MQ, and a = 0.75. 
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Figure 6. Same as Figure 151 but for LBGzSs assuming Afmin 
1.5 X 10"/i-iM(7), Ml = 5 X 10i2/i-iM(7), and a = 0.75. 



two-parameter planes after marginalizing over the remain- 
ing one parameter. The two-dimensional likelihood contours 
represent Ax^ ~ 2.3, 6.17 and 11.8 which, if each bin of 
the correlation function is independent, should correspond 
to 68.3, 95.4 and 99.7% confidence levels (Press et al. 1986). 
Strictly speaking, however, the sampled correlation function 
bins are not completely independent, and these confidence 
levels should be regarded simply as approximate estimates. 

Examine first the parameters for LBGz4s which are 
fairly strongly constrained by the observations (Fig.|3J. Top- 
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Figure 7. The likelihood functions (Ax^) for M^in (dashed lines) 
and for Mi (solid lines) obtained after marginaUzing over the 
other two parameters. Upper and lower panel is for LBGz4 and 
LBGz5, respectively. 



left panel shows the likelihood map on Mmin-Mi plane. 
As pointed out earlier by Berlind & Weinberg (2002), Bul- 
lock et al. (2002), and Moustakas & Somerville (2002), Mi 
and Mniin are mainly constrained by the number density 
and their clustering amplitude on large scales {6 > 1'), re- 
spectively. As Figure|21shows, the observational uncertainty 
in the clustering amplitude for LBGz4s is fairly small and 
5n/n ~ 12%. Thus we have relatively tight constraints on 
those two parameters. The constraint on a is weak because 
of the degeneracy with the other two parameters. However, 
it is clear that the data favor a < 1 implying that the galaxy 
formation is less efficient (or small galaxies merge more effi- 
ciently to form larger ones) in more massive halos. 

Turn next to LBGz5 (Fig. 2J ■ The constraints on this 
population are not so tight because of much larger uncer- 
tainties in the clustering amplitude (Fig. |5J and in the num- 
ber density, 5n/n ~ 62%. Nevertheless the constraints on 
the parameters for LBGzSs seem very similar to those for 
LBGz4, and we are not able to detect any significant differ- 
ence of the parameter values of LBGs between z — 4 and 
5. 

Figures |S] and |S] compare the observed angular two- 
point correlations with the halo model predictions based on 
preferred parameters for LBGz4s and LBGz5, respectively. 
In plotting the model predictions, we adopt Mmin = 1.6 x 
10"/i"^Mq, Ml = 8xlO^^ft"^M0, and a = 0.75 for LBGz4, 
and Mmin = 1.5 x lO^/i'^M©, Mi = 5 x lO^^ft-^M©, and 
a = 0.75 for LBGz5. Given the approximate and empirical 
nature of the halo model, the overall agreement is satisfac- 
tory. 



Table 2. The galaxy-number weighted average mass of hosting 
halo (Miiost) (in units of h~^MQ) and the expected number of 
galaxies per one halo (Ng) for typical values of HOF parameters. 



Sample (Mmin, Mi, a) 


(Ng) 


(Mhalo) 


LBGz4 (1.6 X 10^^ 2.4 x 10^^, 0.5) 
LBGz5 (1.4 x 10", 1.4 x lO^^, 0.5) 


0.38 
0.45 


6.3 X 10" 
4.5 X 10" 



4.2 Characteristics of LBG hosting halos 

Let us look into more carefully the hosting halo masses for 
the two LBG samples. Figure|7|plots the likelihood functions 
(Ax^) for Mmin (dashed lines) and for Mi (solid lines) after 
marginalizing over the remaining two parameters. Clearly, 
both likelihood functions look similar, indicating little evo- 
lution of properties of the hosting halos from z ~ 5 to z ~ 4. 

In order to estimate the characteristic mass of hosting 
halos and the typical number of galaxies per halo, we intro- 
duce the following two quantities: the average mass of the 
hosting halo (member galaxy number weighted): 



(Mhost) 



dM M Afg(M)ni,aio(M) 



and the average number of galaxies per halo: 
Ng{M)ni,MM,z)dM 



. nhalo(M,2)dM 



(14) 



(15) 



Those are evaluated assuming typical sets of HOF pa- 
rameters that we found in the previous subsection (Ta- 
ble I^J. The average mass of the hosting halos for LBGs is 
(5-6)xlO"/i'^M0, and the average number of galaxies per 
halo is ~ 0.4. Thus LBGs have an approximate one-to-one 
correspondence to relatively less massive halos. 



4.3 Evolution of properties of the hosting halos 
for LBGs 

Turn next to the evolution of the hosting halos for LBGs. 
Our current analysis indicates that the minimum halo mass 
Mmin for LBGs is almost the same ~ 1.5 x 10^^ Mq at 
z ~ 4 and z 5. The halo mass accommodating more than 
one galaxy. Mi , seems increasing as time although barely at 
a 1-a level. 

The almost identical analyses by Moustakas & 
Somerville (2001) and Bullock et al. (2002) for 2 ~ 3 
LBGs (Steidel et al. 1998; Adelberger et al. 1998; Adel- 
berger 2000) indicate that Mnin = 1.3 x 10^°/i"^Mq and 
Ml = 6 X 10^^ /i"^M0 for a = 0.8, and that Mmin = 
(0.4 - 8) X lO^^/i-^M©, Ml = (6 - 10) x lO^^ft-^M© and 
0.9 < a < 1.1, respectively. Taking account of the rela- 
tively large uncertainties in those estimates, their results 
are consistent with ours, and indeed the combined results 
may indicate an evolutionary trend of decreasing Afmin and 
increasing Mi with decreasing z. The different selection cri- 
teria at different redshifts may induce an artificial system- 
atic effect in estimating hosting halo mass, but this is not 
the case here. If the limiting flux of a sample is brighter, 
galaxies in the sample have a smaller number density and 
usually a higher clustering amplitude. This leads to increas- 
ing Mmin and Ml simultaneously. However, the limiting ab- 
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solute magnitudes for LBGz3s, LBGz4s and LBGzSs are 
Ml 700 = -19.3 + 5 log ft, Ml 700 = -19.0 + 51og/i and 
Mi7oo ~ —19.7 + 51og/i, respectively. Therefore, it is very 
unlikely that the difference in the limiting magnitude solely 
accounts for the systematic (although weak) trends in Afmin 
and Ml . 

In summary, the hosting halos for LBGs are character- 
ized as follows; (i) Afmin is about ~ 1.5 x lQ"h~^ Mq both 
at z ~ 4 and 2^5, while it decreases to about Mmin = 
(0.4-8) X 10^° ft" ^Mq at 2 ~ 3. (ii) Mi increases with time. 
Ml ~ 1.4 X 10^^ 2.4 X 10^^ and (6 - 10) x 10^^/i"^Mq for 
2 = 5, 4 and 3, respectively. 



5 RESULTS ON LYMAN-ALPHA EMITTERS 

Let us turn to LAEz5s. As shown in Figure m their angu- 
lar correlation function exhibits a somewhat irregular shape. 
This is more clearly seen in the plot of the bias (lower panel 
of Figure |5J. On scales less than 120 arcsec, the bias in- 
creases with decreasing separation similarly to LBGs, and 
its amplitude is in the range between those for LBGz4s and 
LBGzSs. On the other hand, on larger scales the bias factor 
is rather high, which is in a marked contrast with LBGs. The 
number density of LAEs is higher than that of LBGz5s but 
it has a large uncertainty, 5n/n ^ 65%. It should be noted 
that the survey volume of LAEz5 is small, (30/i~^Mpc)^ (co- 
moving volume) , thus it is possible that these measurements 
are significantly affected by the cosmic variance. 

For reference, we give here some numbers which are 
useful in the following discussion; at the redshift of LAEz5s, 
2 ~ 4.86, 1 arcmin corresponds to 1.56/i~^Mpc (comoving), 
and the average number of halos with the mass larger than 
M within the survey volume computed from the halo mass 




Figure 9. Contours show the halo model prediction for the galaxy 
number weighted large scale bias factor defined by eq. |I5V . 



function is iV(> M) = 90, 10 and 1, for M = 3 x 10^\ 
1 x 10^^, and 3 x 10^^/i"^Mq, respectively. 

We apply the same likelihood analysis as performed for 
LBGs in the last section, but we find that our simple HOF 
prescription fails to reproduce simultaneously the observed 
angular correlation function and number density. Indeed, no 
parameter set is found to give a reasonably small . The 
most serious discrepancy is the very high correlation am- 
plitude on scales larger than 120 arcsec. In what follows, 
we present two illustrative examples of failed models, which 
would help to search for a possible solution of the prob- 
lem. For this purpose, we plot halo model predictions for 
the number density of LAEz5s in Figure |S1 in which the 
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Figure 10. Comparison of the observed angular correlation func- 
tion of LAEz5 with the model prediction assuming a = 1.0, 
Mmin = 3.0 X and Mi = 1 X lO^/i-lM© for the 

solid line, and Mi = 1 X lO^^/i"^*/© for the dashed line. The 
dotted line shows the 2-halo term only which does not depends 
on Ml. 
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Figure 11. Same as Figure [TT71 but for different HOF parameter 
sets Mmin = 3.0 X lO^i/i-lMo, a = 1.0, and Mi = 1.0 X lO^^ 
(the solid line) and 10^'' (dashed line). 



gray region indicates the 1-cr range of the observed number 
density, and halo model predictions for the large scale bias 
factor defined by eq. © in Figure |U] 

The first example is as follows; taking the amplitude of 
the correlation function on large scales (6 = 7 — 9 for 9 > 120 
arcsec, see Figure I^J, one finds in Figure |^ that Mmin = 



(2 — 5) X W^^ Mq is required irrespective of a. This value 
combined with the observed number density of LAEz5s gives 
roughly Mi ~ 10^^h~^MQ (or less) for a reasonable range 
of a (Fig.lHJ. This parameter set, however, predicts a much 
higher correlation amplitude than observed on smaller scales 
(Fig. Iim . If one attempts to have a reasonable correlation 
amplitude on smaller scale, a very large Mi (more than 
W^'^ Mq at least) is required (see Fig. 1101 . which leads 
to too small a number density on the other hand. 

One possible way to reconcile this discrepancy is to 
modify the halo model that we have adopted. Figure^] in- 
dicates that it is the 1-halo term that boosts the correlation 
amplitude on smaller scales. Therefore, our assumption in 
the halo model that the galaxy distribution follows the dark 
matter distribution may not hold for LAEzSs. In this pa- 
per, however, we do not attempt to develop the HOF model 
by allowing a possible variation on the galaxy distribution 
within halos, because the statistical accuracies of the correla- 
tion functions on small scales are still relatively low. Indeed, 
the number of small separation pairs is very small; the num- 
ber of LAE pairs which fall into the smallest separation bin 
is two, and that into the second bin is five. This means that 
halos having more than one LAE are limited, and one has 
to keep in mind that the correlation function measurements 
on these scales are based on LAEs in such a small number 
of halos. Given a possible large uncertainty in the measure- 
ments of small-scale clustering, it is premature to explore 
an alternative model in great detail at this point. We also 
note that given the limited number of LAEzSs, it is not to 
be denied that their correlation signals are contaminated by 
the presence of the possible foreground galaxies in one or a 
few clusters. We must wait for future data extending over a 
much larger survey area. 

In addition to the above problem, there is another prob- 
lem concerning the predicted number of halos in the survey 
volume. If one takes Mmin = 3 x 10^^ Mq, the expected 
mean number of halos with mass larger than this value in the 
survey volume is iV(> 3 x 10^^/i"^Mg) ~ 1. This suggests 
that there are only a few hosting halos in the survey volume. 
If this picture is right, it is unlikely that the 2-halo term is 
correctly measured from such a small number of halos. 

The above consideration implies that a simple halo 
model may not be readily applicable to the LAEs. Note 
that the parameter set {Mmin, Mi, a) = (3 x 10^^/i~"^Mq, 
1 X lO"ft-^M0, 1.0) gives (Ng) = 45 and (Mhaio) = 
5.2 X 10^^ Mq. Therefore, if real, we can conclude at least 
that the nature of hosting halos and the relation between 
galaxies and halos for LAEs are very different from those 
for LBGs. 

The second example is as follows. Let us ignore the data 
on large scales {9 > 120 arcsec) for a while, and adopt the 
value of b around 1 arcmin, i.e., fe ~ 4. This gives Mmin ~ 
(2-4) X lO"/i"^M0 (Fig.EJ. This value combined with the 
observed number density gives roughly Mi ~ lO^^fe-^M© 
for a reasonable range of a (Fig.|Sll. This parameter set pre- 
dicts a slightly higher correlation amplitude on small scales 
than observed as shown in Figure lTTl However, this discrep- 
ancy may not be taken so seriously, because of the limited 
statistical significance as mentioned above. Better agree- 
ment is obtained by setting larger Mi ^ lO^"^ Mq (Fig. 
IllH . which, however, predicts a smaller number density than 
observed. The discrepancy in the number density becomes 
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smaller if one takes a smaller a (Fig. [HJ. If we take {Mmin, 
Mi)=(3 X 10"/i"^Mq, 1 X W^^h'^MQ), the characteristic 
values are (Ng) = 0.67 and (Mhaio) = 9.0 x lO^/i'^M© for 
a = 0.75, and {Ng) = 0.75 and (Mhaio) = 7.8 x IO^/i-^Mq 
for a = 0.5. These values are similar to those for LBGs. 

The above model provides acceptable agreement with 
both the correlation function on scales smaller than 120 arc- 
sec and the number density. The predicted correlation func- 
tion on larger scales, however, has a much lower amplitude 
than observed. A possible interpretation to this discrepancy 
is that the current survey volume is one of the overdense re- 
gions on large scales by chance, and the LAEs in the region 
have accidentally acquired the high correlation amplitude. 
In fact, Shimasaku et al. (2003) have extended the survey 
area of LAEz5s to the north and found a high overdensity 
of LAEz5s over a circular region of 5 arcmin (8/i~^Mpc) ra- 
dius. They have suggested that it may be a progenitor of 
a present-day massive cluster of galaxies. Half of this circu- 
lar region is inside our survey area. LAEz5s associated with 
this large-scale overdense region will have an unusually high 
correlation amplitude on large scales. If this is the case, the 
HOF parameter values obtained ignoring the large-scale cor- 
relation function will be close to the typical values of LAEs 
at z ~ 5. 



6 SUMMARY AND DISCUSSIONS 

We have analyzed three high-redshift galaxy samples created 
from the Subaru Deep Field (SDF) survey data; LBGs at 
z ~ 4 (LBGz4s), LBGs at z ~ 5 (LBGz5s) and LAEs at 
z ~ 4.86 (LAEz5s), and explored the implications of their 
number density and angular clustering in the framework of 
the halo occupation function (HOF). 

Our major findings are summarized as follows; 

(i) The two LBG samples can be well described by the 
halo model with an appropriate HOF in an approximate 
fashion. 

(ii) The hosting halos for LBGz4s and LBGz5s are more 
massive than Mmin ~ 1.5 x lO"/i"^M0. Since the expected 
number of LBGs per halo with M > Mmin is ~ 0.5, there 
is an approximate one-to-one correspondence between halos 
and LBGs. This is basically consistent with the results pre- 
viously found for LBGs at z ^ 3 (Mo, Fukugita 1996; Steidel 
et al. 1998; Jing & Suto 1998 Moustakas & Somerville 2001; 
Berlind & Weinberg 2002; Bullock et al. 2002). 

(iii) On the other hand, this may also indicate that a 
large fraction of dark halos do not host LBGs brighter than 
Ml 700 ~ — 19 mag. Nevertheless this does not necessarily 
mean that there is no galaxy in such halos. Franx et al. 
(2003) have found a population of red galaxies at 2 ~ 3 
which do not have active star formation and so may not be 
easily detectable by the Lyman break technique because of 
the very faint UV continuum emission. They estimate that 
the number density of such red galaxies is about half that of 
LBGs at the same redshift. Thus it may be the case that a 
fraction of the dark halos at z ~ 4 — 5 host such red galaxies 
rather than bright LBGs that we have discussed here. 

(iv) The LBG samples at z ~ 3, 4 and 5 discussed here 
have very similar limiting absolute magnitudes, M1700 ~ 
— 19. On the other hand, the minimum mass of their hosting 
halos seems decreasing with time, although its statistical 



significance is not strong. If true, this means that the star 
formation efficiency per unit dark matter mass, L1700 /Mhaio, 
increases with time. This increase may suggest that cold gas 
gradually accumulates in LBGs, if the star formation rate is 
simply proportional to the amount of cold gas. 

(v) There is a weak indication for Mi to increase slightly 
with time in the LBG samples. If this is indeed the case, it 
may be explained by the mutual merging of LBGs. 

(vi) For LAEz5s, our simple HOF prescription fails to 
reproduce simultaneously the observed angular correlation 
function and number density. No parameter set gives a rea- 
sonably small ■ This is mainly because the LAEz5s exhibit 
very strong clustering signal on scales larger than 120 arcsec. 

In fact, the nature of LAEz5s is still uncertain in the 
current result; models which match both the high correla- 
tion amplitude on large scales and the number density of 
LAEz5s predict much higher correlation amplitude on small 
scales than observed. A possible interpretation of this dis- 
crepancy is that the distribution of LAEs within halos differs 
from that of dark matter. If this is the case, the simple halo 
model description for the LAE we adopted in this paper 
needs to be improved. Also the observational indication for 
the discrepancy should be carefully confirmed with future 
larger and more accurate data samples. 

Alternatively, if one constructs models which match 
both the correlation function on small scales and the number 
density of LAEz5s, they imply a lower correlation amplitude 
on large scales than observed. This may be simply due to 
statistical fluctuation in a sense that the current data do 
not represent a fair sample of the universe at z ~ 5 as indi- 
cated by by Shimasaku et al. (2003). The HOF parameters 
derived from the fit to the observed data except for the large 
scale correlation function are M„,„ ~ 3 X lO"ft"^Af0 and 
All ~ 1 X lO^^/i~^M0 for a reasonable range of a. This gives 
(Ng) ~ 0.7 and (Mhaio) ~ 8 x lO^^/i-^M© for a = 0.75, sug- 
gesting an approximate one-to-one correspondence between 
LAEs and halos as found for LBGs. 

In order to distinguish the above two pictures, we need 
a much larger observational sample after all. 



ACKNOWLEDGMENTS 

T.H. and M.O. acknowledge support from Japan Society 
for Promotion of Science (JSPS) Research Fellowships. I. K. 
gratefully acknowledges support from the Takenaka-Ikueikai 
fellowship. This research was also supported in part by the 
Grants-in-Aid from Monbu-Kagakusho and Japan Society 
of Promotion of Science (12640231, 13740150, 14102004, and 
1470157). Numerical computations presented in this paper 
were carried out at ABAC (the Astronomical Data Analysis 
Center) of the National Astronomical Observatory, Japan 
(project ID: mys02a, yys08a). 



REFERENCES 

Adelberger K. L., 2000, in ASP Conf. Ser. 200, Clustering 
at High Redshift, ed. A. Mazure, O. Le Fevre, & V. Le 
Brun (San Francisco: ASP), 13 

Adelberger K. L., Steidel C. C., Giavalisco M., Dickinson 
M., Pettini M., Kellogg M., 1998, ApJ, 505, 18 



© 2003 RAS, MNRAS 000.1110 



Properties of host haloes of Lyman-break galaxies and Lyman- alpha Emitters 11 



Bardeen J. M., Bond J. R., Kaiser N., Szalay A. S. 1986, 
ApJ, 304, 15 

Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 291 
Benson A. J., Cole S., Frenk C. S., Baugh C. M., Lacey 
C. G., 2000, MNRAS, 311, 793 
Berlind A. A., Weinberg D. H., 2002, ApJ, 575, 587 
Bond J. R., Cole S., Efstathiou C, Kaiser N., 1991, ApJ, 
379, 440 

Bower R. C, 1991, MNRAS, 248, 332 
Bruzual A. C, Chariot, S., 1993, ApJ, 405, 538 
Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., 
Kravtsov A. V., Klypin A. A., Primack J. R., Dekel A. 
2001, MNRAS, 321, 559 

Bullock J. S., Wechsler, R. H., Somerville R. S., 2002, MN- 
RAS, 329, 346 

Cooray A., Hu W., Miralda-Escude J. 2000, ApJ, 535, L9 

Cooray A., Sheth R., 2002, Phys. Rep., 372, 1 

Cowie L. L., Hu E. M., 1998, AJ, 115, 1319 

Furusawa, H., Shimasaku, K., Doi, M., Okamura, S., 2000, 

ApJ, 534, 624 

Franx, M. et al., 2003, ApJ, 587, L79 
Giavalisco M., Dickinson M., 2001, ApJ, 550, 177 
Hamana T., Yoshida N., Suto Y., 2002, ApJ, 568, 455 
Hamana T., Yoshida N., Suto Y., Evrard A. E., 2001, ApJ, 
561, L143 

Hu E. M., Cowie L. L., McMahon R. C, 1998, ApJ, 502, 
L99 

Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., 

Cole S., Evrard A. E., Couchman H. M. P., Yoshida N., 

2001, MNRAS, 321, 372 

Jing Y. P., 1998, ApJ, 503, L9 

Jing Y. P., Mo H. J. Boner C, 1999, ApJ, 494, 1 

Jing Y. P., Suto Y., 1998, ApJ, 494, L5 

Kauffmann G., Colberg J. M., Diaferio A., White, S. D. M., 

1999, MNRAS, 303, 188 

Komatsu E., Kitayama T. 1999, ApJ., 526, LI 

Lacey C, Cole S. 1993, MNRAS, 262, 627 

Limber D. N. 1953, ApJ, 117, 134 

Ma C.-R, Fry J. N. 2000, ApJ, 543, 503 

Madau P., Ferguson H. C, Dickinson M. E., Giavalisco M, 

Steidel C. C, Fruchter A., 1996, MNRAS, 283, 1388 

Malhotra, S. & Rhoads, J. E. 2002, ApJ, 565, L71 

McClelland J., Silk J. 1977, ApJ, 217, 331 

Mo, H. J., Fukugita, M. 1996, ApJ, 467, L9 

Mo, H. J., Mao, S, White, S. D. M. 1999, MNRAS, 304, 

175 

Mo, H. J., White, S. D. M. 1996, MNRAS, 282, 347 
Moustakas L. A., Somerville R. S., 2002, ApJ, 577, 1 
Navarro J., Frenk C, White S. D. M. 1996, ApJ, 462, 563 
Navarro J., Frenk C, White S. D. M. 1997, ApJ, 490, 493 
Neyman J., Scott E. L. 1952, ApJ, 116, 144 
Ouchi, M. 2003, Ph.D. Thesis, University of Tokyo 
Ouchi M., et al., 2001, ApJ, 558, L83 
Ouchi M., et al., 2003a, ApJ, 582, 60 

Ouchi M., et al., 2003b, Carnegie Observatories 
Astrophysics Series, Vol. 3: Clusters of Galax- 
ies: Probes of Cosmological Structure and Galaxy 
Evolution , ed. J. S. Mulchaey, A. Dressier, and 
A. Oemler (Pasadena: Carnegie Observatories, 
|http : / /www.ociw.edu /ociw /sy mposia /ser ies/symposium3 
/proceedings. html) 

Ouchi M., et al., 2003c, ApJ submitted tastro-ph/0309655| l 



Ouchi M., et al., 2003d, ApJ submitted | |astro-ph/0309657t 
Peacock J. A., Dodds S. J., 1996, MNRAS, 280, L19 
Peacock J. A., & Smith R. E. 2000, MNRAS, 318, 1144 
Peebles P. J. E. 1974, A&A, 32, 197 

Peebles P. J. E., 1980, The Large-Scale Structure of the 
Universe Princeton Univ. Press, Princeton, NJ 
Porciani C, Giavalisco M., 2002, ApJ, 565, 24 
Press W. H., Flannery B. P., Teukolsky S. A., Vetterling 
W. T., 1986, Numerical Recipes (Cambridge: Cambridge 
Univ. Press) 

Press W. H., Schechter P., 1974, ApJ, 187, 425 
Schaerer D., 2003, A&A, 397, 527 

Scoccimarro R., Sheth R. K., Hui L., Jain B., 2001, ApJ, 
546, 20 

Seljak U. 2000, MNRAS, 318, 203 

Sheth R. K., Jain B.. 1997, MNRAS, 285, 231 

Sheth R. K., Mo. H. J., Tormen G., 2001, MNRAS, 323, 1 

Sheth R. K., Tormen G., 1999, MNRAS, 308, 119 

Sheth R. K., Tormen G., 2002, MNRAS, 329, 61 

Shimasaku K. et al., 2003, ApJ, 586, Llll 

Shimizu, M., Kitayama, T., Sasaki, S., Suto, Y., 2003, ApJ, 

590, 197 

Shu, C, Mao, S., Mo, H. J., 2001, MNRAS, 327, 895 

Steidel C. C, Giavalisco M., Pettini M., Dickinson M., 

Adelberger K. L., 1996, ApJ, 461, L17 

Steidel C. C, Adelberger K. L., Dickinson M., Giavalisco 

M., Pettini M., Kellogg M., 1998, ApJ, 492, 428 

Steidel C. C, Adelberger K. L., Giavalisco M., Dickinson 

M., Pettini M., 1999, ApJ, 519, 1 

Takada M., Hamana T., 2003, MNRAS in press 
(astro-ph/0305381) 

Takada M., Jain B., 2003, MNRAS, 340, 580 
White M., Hernquist L., Springel V., 2002 ApJ, 580, 634 
Yoshikawa K., Taruya A., Jing Y.-R, Suto Y., 2001, ApJ, 
558, 520 



APPENDIX A: DEPENDENCES OF HOF 
PARAMETERS ON THE SHAPE OF THE 
TWO-POINT CORRELATION FUNCTION 

Here we summarize major characteristics of changes in the 
shape of the two-point correlation functions made by vary- 
ing HOF parameters. Figure IA1IA3I show the angular cor- 
relation functions for variety of HOF models. Each Figure 
demonstrates the effect of varying one HOF parameter with 
fixing other two parameters. In all plots, the selection func- 
tion of LBGz4 is adopted for an illustrative purpose. It 
should be noticed that a degree of change in the shape of 
the correlation function made by varying one parameter de- 
pends on a choice of other two parameters, thus the plots 
should be understood as illustrative examples, and we just 
focus on major characteristics from a qualitative point of 
view. 

Fieure lXTI shows the effect of varying A/i on uj. Since the 
2-halo term does not depend on Mi , the correlation function 
on larger scales is not affected by the change of Afi. The 
amplitude of the 1-halo term decreases with increasing Mi . 
This is mainly due to the decrease in the contribution from 
smaller mass halos. 

Figure IX2I is for effect of varying Mmin- The amplitude 
of the 2-halo term increases with Mmin, because of a larger 
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Figure Al. Dependence of varying M\ on the shape of angular 
two-point correlation functions. HOF parameters for three cases 
are inserted in the plot. Two thin lines show 1-halo (dominates 
on small scales) and 2-halo term (dominates on large scales) and 
thin lines show the sum of them. We take the selection function of 
LBGz4 to compute those correlation functions for an illustrative 
purpose. 
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Figure A3. Same as Figure IXn but for a. 



halo term. The change in the slope of the 1-halo term with 
varying a is explained as follows: Since, roughly speaking, 
galaxies in larger halos contribute to the 1-halo term on rel- 
atively larger scales, while galaxies in smaller halos can only 
contribute to the 1-halo term on smaller scales. Varying a 
changes the fraction in the contribution from larger halos to 
smaller mass halos, and accordingly changes the slope of the 
1-halo term. A smaller a results a steeper slope, because of 
a more weight for galaxies in smaller halos. 
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Figure A2. Same as Figure IXTl but for Mmin 



bias factor for more massive halos. The amplitude of the 
1-halo term also increases with Mmin- This is largely the 
decrease in the contribution from smaller mass halos. 

Finally Figure IXTl shows the effect of varying a on u. 
Varying a changes the fraction of galaxies in massive halo to 
less massive halos. A larger a gives more weight for galax- 
ies in massive halos. Since a stronger bias factor for more 
massive halo, a larger a gives a larger amplitude of the 2- 
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